Remanence of Ni nanowire arrays: influence of size and labyrinth magnetic structure 
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The influence of the macroscopic size of the Ni nanowire array system on their remanence state 
has been investigated. A simple magnetic phenomenological model has been developed to obtain 
the remanence as a function of the magnetostatic interactions in the array. We observe that, due to 
the long range of the dipolar interactions between the wires, the size of the sample strongly influence 
the remanence of the array. On the other hand, the magnetic state of nanowires has been studied by 
variable field magnetic force microscopy for different remanent states. The distribution of nanowires 
with the magnetization in up or down directions and the subsequent remanent magnetization has 
been deduced from the magnetic images. The existence of two short-range magnetic orderings with 
similar energies can explain the typical labyrinth pattern observed in magnetic force microscopy 
images of the nanowire arrays. 
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I. INTRODUCTION 

During the last decade, regular arrays of magnetic 
nanoparticles have been deeply investigated^^ Besides 
the basic scientific interest in the magnetic properties of 
these systems, there is evidence that they might be used 
in the production of new magnetic devices and partic- 
ularly in recording mediai^i^ Different geometries have 
been considered, including dots, rings, tubes and wires. 
Recent studies on such structures have been carried out 
with the aim of determining the stable magnetized state 
as a function of the geometry of the particles/^i^ In par- 
ticular, the study of highly ordered arrays of magnetic 
wires with diameters typically in the range of tens to hun- 
dreds of nanometers is a topic of growing interest i^i^iiSiii^ 
Anodization processes to achieve self-ordered nanopores 
in membranes have been proven to be a direct, sim- 
ple, and nonexpensive technique in fabricating templates 
for highly ordered densely packed arrays of magnetic 
nanowires. "'^^ The high ordering, together with the mag- 
netic nature of the wires, gives rise to outstanding co- 
operative properties of fundamental and technological 
interest, since it can determine the success of patterned 
media in high-density information storage. 

Effects of interparticle interactions are, in general, 
complicated by the fact that the dipolar fields depend 
on the magnetization state of each element, which, in 
turn, depends on the fields due to adjacent elements. 
Therefore, the modeling of interacting arrays of wires is 
often subject to strong simplifications such as, for ex- 
ample, modeling the wire using a one-dimensional classi- 
cal Ising modelJ^ii^ Also micromagnetic calculation o^^i^^ 
and Monte Carlo simulationsii have been developed. 
However, these methods typically permit us to consider 
only an array of a reduced number of nanowires, a sit- 
uation far from the state of a regular array, as stated, 
for example, by Sampaio et al.,^^ who observed modifi- 



cations of the remanent magnetization as a function of 
the number of wires in systems with 2-500 elements. 

The magnetization of ferromagnetic nanowire arrays 
has already been studied by magnetic force microscopy 
(MFM) that, in addition, enables us to gain direct mag- 
netic information of individual nano-objects. In previous 
works, MFM measurements have been carried out by ap- 
plying magnetic fields on magnetized and demagnetized 
samples both to study the switching behavior of indi- 
vidual nanowires and to obtain the hysteresis loops of 
the nanowire arrays i^i^iiS In the equilibrium state, the 
nanowires exhibit a homogeneous magnetization along 
the axial direction (with the magnetization of each wire 
pointing up or down). Then, it appears reasonable to 
investigate the role of interactions in a microscopic ar- 
ray using a model of single-domain structures including 
length corrections due to the shape anisotropy. In this 
paper we perform an experimental and theoretical study 
to understand the role of the size of the system on the 
remanence of a Ni nanowire array. We also investigate 
the pattern domain structure of an array, which can be 
explained by considering dipolar interactions in a typical 
hexagonal cell. 



II. EXPERIMENTAL DETAILS 

Self-assembled nanopores with hexagonal symmetry 
have been obtained in alumina matrix by two-step an- 
odization process. Subsequently, pores are filled with 
Ni by electroplating process. Full details of prepara- 
tion method can be found elsewhere. Ni nanowires 
{d — 2R ~ 180 nm in diameter and L = 3.6 /im in length) 
are arranged in a hexagonal pattern with D = 500 nm 
lattice constant. 

The hysteresis loops along the axial direction has been 
measured by a superconducting quantum interference de- 
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vice (SQUID). The local magnetization distribution has 
been studied using a MFM equipment from Nanotec 
Electronica"^*^. Such a system, working in noncontact 
mode, allows us to acquire simultaneously the topogra- 
phy of the surface and the magnetic force gradient map. 
The MFM system has been conveniently modified in or- 
der to apply continuously a magnetic field of ±0.2 T along 
the in-plane direction and pulses along the out-of-plane 
directionJ^ By using this so-called variable field mag- 
netic force microscopy technique, the reversal process of 
Ni nanowires has been studied. In order to avoid the 
tip influence on the magnetic state of the nanowires^ 
MESP low moment MFM probes have been used. 



of one of the wires. In fact, Hani can be recognized as 
the value of the coercivity He of each individual wire and 
can be calculated as 

Ham = H,^ aMo [N, - N,) (3) 

with < a < 1 a factor determined by the magneti- 
zation reversal mechanism and dipolar interactions be- 
tween the wires in the arrayi^. The demagnetizing fac- 
tors are given bySe = (1/2) F21 [AR'^/L^] - {AR/SnL) 
and N, = l^ F21 [AR^/L^] + {8R/3ttL), where F21 [x] = 
F21 [—1/2, 1/2, 2, —x] is a hypergeometric function. 



III. MONTE CARLO SIMULATIONS 

As a consequence of the large aspect ratio of the wires 
investigated, L/R = 40, the anisotropy they present is 
mainly shape anisotropy. In this case, the individual 
wires can be considered as nearly single-domain struc- 
tures with two stable states: the magnetic moment point- 
ing up or down. However, the behavior of the ar- 
ray as a whole differs from a pure bistable magnetic 
state due to the magnetostatic interactions between the 



nanowires 
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In order to model the hysteresis loop of 
the array, we develop Monte Carlo simulations consider- 
ing magnetostatic interactions. The starting point of the 
model assumes that each nanowire has a magnetization 
oriented along any of the two axial directions (z axis) due 
to the shape anisotropy and all the wires in the array in- 
teract magnetostatically. The internal energy, E, of the 
array with N identical wires can be written as 



E = ^ioMoV 




DijaiUj - {Ha 



N 

Hani) ^ Cr.j 
i=l 



(1) 

where Mq is the saturation magnetization and V — ttR^L 
is the volume of each wire. The variable Ui can take the 
values ±1 on a site z of a two-dimensional array, allowing 
the magnetic nanowires to point up {ui = -1-1) or down 
(cTj = — 1) along the axis of each individual wire. The first 
term in the above equation is the dipolar interaction of 
all pairs of magnetic wires. The coupling constant Dij 
has been calculated by Laroze et al.^^ and is given by 
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with Tij the distance between the magnetic wires at sites 
i and j. Note that since Dij is positive, the dipolar in- 
teraction favors an antiparallel alignment between the 
magnetic nanowires. The second term in Eq. [1] corre- 
sponds to the contribution of an external magnetic field. 
Ha, applied along the axis of the wire and the third term, 
Hani, corresponds to the field representing the magnetic 
shape anisotropy of a single wire, i.e., the reversal field 



IV. RESULTS AND DISCUSSION 
A. Hysteresis loops 

Figure 1(a) shows the hysteresis loop of the Ni 
nanowire array along the axial direction measured with 
SQUID magnetometer. The size of the sample is 16 mm^. 
The magnetic information obtained from the hysteresis 
loop suggests a system with well defined easy axis par- 
allel to the nanowires due to the shape anisotropy. The 
coercive field is He = 215 Oe and the remanent magne- 
tization M* = 0.379 Mo, with Mq = 480 emu/cm^. For 
the wires under consideration, Nx — = 0.469. 

In order to understand the effect of the dipolar inter- 
actions on the hysteresis loop, we performed numerical 
simulations for the same array. From the measured value 
for He, we obtain Hani = 2850a Oe through Eq. (3) 
which defines a « 0.08. By considering this value into the 
energy expression [Eq. (1)], we are in conditions to sim- 
ulate the hysteresis loop. Monte Carlo simulations were 
carried out using Metropolis algorithm with local dynam- 
ics and single-spin-flip methods.— The initial state at Ha 
— 2.0 kOe, higher than the saturation field, considers the 
magnetization of all the wires aligned with the external 
field. The field was then linearly decreased at a rate of 
300 Monte Carlo steps for AH = 0.01 kOe. The new ori- 
entation of the magnetic nanowire was chosen arbitrarily 
with a probability p = min [1, ea;p(— AiJ/fesT)], where 
AE is the change in energy due to the reorientation of 
the wire, and ks is the Boltzmann constant. Figure 1(b) 
illustrates the hysteresis loop for samples with size rang- 
ing from 4 fim^ {N = 14) to 2116 fim'^ {N = 9699). Each 
loop is the result of the average of five independent real- 
izations. Because in our calculations the internal struc- 
ture of the wire is not considered, the coercivity has a fix 
value independent of the number of wires in the array. 
From this figure, we can observe that the size of the sam- 
ple strongly influences the shape of the loop as a whole 
and the particular role of the remanence. With present 
standard computational capabilities, it is not possible to 
obtain hysteresis loops with N higher than 10000. To 
describe the remanence of bigger arrays such as the ones 
experimentally investigated, we propose an alternative 
approach presented in next section. 
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FIG. 1: (a) Hysteresis loops as measured by SQUID with 
the external field applied parallel to the wire direction. The 
gray dots correspond to the magnetic moment deduced from 
the MFM images (see Sec. IV C). (b) Hysteresis loops per- 
formed by numerical simulations considering different sizes of 
the sample. 



B. Remanent magnetization 

In order to understand the role of the size of 
the sample, we calculate the dipolar interaction en- 
ergy per wire for arrays of different sizes, Eint = 
Et7'Ef=.+i (A,/MoiV). By inspecting Fig. 1(b), we 
observed that it is possible to introduce a phenomenolog- 
ical analytical function that allows us to obtain the re- 
manent magnetization as a function of the magnetostatic 
interaction present in the array. The reduced remanence 
rur can be written as 



where < rUr < 1. The dependence of calculated rema- 
nence rrir (dotted line) and interaction energy E^nt (solid 
line) is depicted in Fig. 2 as a function of the size of the 
sample (i.e., number of nanowires N). Black dots illus- 
trate rrir obtained from the hysteresis loops shown in Fig. 
1(b). For an array of 14400 fim^ {N « 70000), the rema- 
nent magnetization converges to the value rn^ = 0.38, 
which is in excellent agreement with the remanence ob- 
tained with the SQUID measurements. These calcula- 
tions establish a lower limit to the number of wires that 
have to be used in simulations to reach the experimental 
remanence value. Nevertheless, experimentally, it is not 
possible to measure arrays less than a few mm^ in surface 
area with the required precision. At present, samples for 
experiments are 16 mm^, a size beyond saturation. Con- 
sequently, the theoretically established dependence of re- 
manence on the sample size of such nanostructures is not 
yet experimentally confirmed. However, in the future, 
if smaller samples can be measured, it is important to 
consider this lower size limit. 
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FIG. 2: Magnetostatic interaction per wire in the array (solid 
line) and remanent magnetization (dotted line) as a function 
of the size of the sample. Black dots correspond to the rema- 
nent magnetization obtained from hysteresis loops simulated 
in Fig. 1(b). 

Now we further investigate the validity of Eq. 2] by 
calculating the remanence for different samples we found 
in the literature. Assuming that the size of the sample 
is big enough to reach convergence of the remanence, we 
use measured values of He to calculate a and rrir through 
Eqs. (3) and (4). In order to assure convergence, we fix 
N — 70000. Table I summarizes the geometrical parame- 
ters of the array. He, a (which accounts for the influence 
of magnetic interactions in the reversal process and coer- 
civity), the measured m*, and calculated values of re- 
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manence. Note the agreement between experimental and 
calculated values through Eq. (4). Deviations between 
SQUID measurements and analytical results originated 
from the dispersion of the lengths and positions of each 
wire in the array and a reduction in the homogeneity of 
the diameter of nanopores.-iS Notice that in the sample 
defined by d = 55 nm and D = 65 nm, the wires are 
very close and then strong interactions are present be- 
tween contiguous nanowires. Due to this interaction, the 
remanence decreases as evidenced in the measured and 
calculated values of the remanence in Table I. 



TABLE I: Parameters for different Ni nanowire arrays. Geo- 
metrical parameters, He, and m* have been measured in this 
paper [superscript 1] and taken from Ref. 28 [superscript 2] 
and Ref. 8 [superscript 3]. 
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0.35 
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C. The patterned domain structure 

MFM images have been obtained in different remanent 
states, since our MFM system permits us to apply mag- 
netic fields in the course of the microscope operation. 

The images in Fig. 3 illustrate different patterns ob- 
tained after applying increasing axial magnetic field. It 
is worth to mention that the MFM tip was previously 
saturated with positive field, in the same direction as 
the subsequent in situ applied magnetic field. The white 
contrast corresponds to nanowires with the magnetiza- 
tion oriented opposite to the tip field direction. When 
the magnetization of the nanowires points in the tip field 
direction, we obtain black contrast. The initial state [Fig. 
3(a)] was achieved after magnetic saturation of the sam- 
ple with a negative magnetic field. The images obtained 
after applying in situ magnetic fields of 90, 190, 290, 400, 
and 500 Oe along the axial direction (see Figs. 3(b)-3(f), 
respectively) show us the evolution of the magnetic state 
of the individual nanowires. Notice the increment of the 
number of black nanowires after applying increasing mag- 
netic field. 

In a previous report, quantitative information of net 
magnetization from MFM images was analyzed. By 
counting the number of wires pointing in each direction. 




(d) Ha = 300 Oe (e) Ha = 400 Oe (f) Ha = 500 Oe 



FIG. 3: MFM images of the (a) initial state and after applying 
fields of (b) 90 Oe, (c) 190 Oe, (d) 290 Oe, (e) 400 Oe, and 
(f) 500 Oe parallel to the tip field. 



the remanent magnetization value, TO*can be obtained as 

. N,„ - Nb 



(5) 



where N^^ and Nf, are the numbers of wires with mag- 
netization pointing up and down, respectively. Counting 
the black and white points in the MFM images in Fig. 3, 
m* was obtained for different values of Ha, which are il- 
lustrated with dots in Fig. 1(a). In the MFM results, the 
effect of the stray field of the tip must be taken into ac- 
count. Moreover, since the images have been acquired in 
remanent states, the calculated magnetization values are 
slightly lower than the data obtained from the SQUID 
measurements. 

Interactions play a fundamental role on the magnet- 
ically patterned structure of the samples. The pat- 
terned structure in an array, in principle, obeys to an 
antiferromagnetic-like alignment due to the magnetic in- 
teraction between the nanowires. As earlier reported^^ 
for a square lattice, each of the four nearest neighbors 
aligns antiparallel and the magnetic structure of the array 
exhibits a checkerboard pattern. However, when we con- 
sider a typical hexagonal cell, as in Ref. 15, we have two 
almost degenerate states. At T = 0, the configuration in 
Fig. 4(a) has, for d — 180 nm, D ~ 500 nm and L ~ 
3.6 /im, a 10 % less energy than the configuration in Fig. 
4(b). Due to such a small difference, the temperature, 
lattice disorder, or the magnetic history of the sample 
allows the array to exhibit any of both short-range con- 
figurations. Then, in a regular array, a mixture between 
both states is observed which originates the labyrinth 
pattern shown in Fig 4(c). This figure has been obtained 
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by means of Monte Carlo simulations, as explained be- 
fore, starting from a saturated sample and decreasing the 
external field until the coercive value. In this state, al- 
most the same number of wires has their magnetization 
pointing up (white in Fig. 4) or down (black in Fig. 4), 
and nearest-neighbor parallel magnetic moments are or- 
ganized in structures such as the ones shown in Fig. 4(a) 
and 4(b). 



a) 0) 




FIG. 4: Magnetic configuration of a typical hexagonal cell 
with (a) minimum energy and (b) first excited state. Black 
(white) dots represent a wire with its magnetization pointing 
up (down). The energy difference between both configura- 
tions is 10 %. (c) Simulated patterned domain structure at 
remanence state. Image size: 15 fj,m x 15 //m. 

A comparison between simulated and MFM labyrinth 
images confirms the at least qualitative agreement of this 
approach. 

V. CONCLUSIONS 

In conclusion, by means of theoretical studies and ex- 
perimental measurements, we have investigated the im- 



portant role of magnetostatic interaction in the magnetic 
properties of nanowire arrays. We have derived an an- 
alytical expression that allows one to obtain the rema- 
nent magnetization as a function of the magnetostatic 
interactions presented in the array. Our results lead us 
to conclude that, because of the long-range order of the 
dipolar interactions between the wires, the size of the 
sample strongly influences the remanence of the array. 
In order to guarantee reproducibility, it is important to 
consider a sample which contains a minimum of 70000 
wires. Also, the typical labyrinth pattern observed in 
the MFM images has been explained by a simple model 
considering the presence of two magnetic patterns of the 
basic cell of an hexagonal array. The MFM proves to 
be a useful method in studying the reversal magnetiza- 
tion process in nanostructures. Moreover, this powerful 
technique allows us to observe the individual evolution of 
the magnetic state of hundreds of nanowires under an ex- 
ternal magnetic field. Good agreement between SQUID 
and MFM measurements and theoretical simulations is 
obtained. 
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